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ABSTRACT 

Using mathematical formalism borrowed from dynamical systems theory, a complete ana- 
lytical investigation of the critical behaviour of the stationary flow configuration for the low 
angular momentum axisymmetric black hole accretion provides valuable insights about the 
nature of the phase trajectories corresponding to the transonic accretion in the steady state, 
without taking recourse to the explicit numerical solution commonly performed in the litera- 
ture to study the multi-transonic black hole accretion disc and related astrophysical phenom- 
ena. Investigation of the accretion flow around a non rotating black hole under the influence 
of various pseudo-Schwarzschild potentials and forming different geometric configurations 
of the flow structure manifests that the general profile of the parameter space divisions de- 
scribing the multi-critical accretion is roughly equivalent for various flow geometries. How- 
ever, a mere variation of the polytropic index of the flow cannot map a critical solution from 
one flow geometry to the another, since the numerical domain of the parameter space re- 
sponsible to produce multi-critical accretion does not undergo a continuous transformation 
in multi-dimensional parameter space. The stationary configuration used to demonstrate the 
aforementioned findings is shown to be stable under linear perturbation for all kind of flow 
geometries, black hole potentials, and the corresponding equations of state used to obtain the 
critical transonic solutions. Finally, the structure of the acoustic metric corresponding to the 
propagation of the linear perturbation studied are discussed for various flow geometries used. 
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1 INTRODUCTION 

Astrophysical blackholes manifest their presence only gravitationally. No spectral information can directly be obtained from these candidates 
because of the presence of the event horizon. One can only, therefore, rely on accretion processes to understand their observational signatures 
l Pringlelll98ll ; (tCato et alJI 19981 ; iFrank et alj|2002h . At large distances from the accretor, black hole accretion is usually subsonic. The inner 



boundary condition imposed by the event horizon is determined by the requirement that the flow will be of a supersonic nature very close to 
the accretor. Black hole accretion, thus, usually demonstrates transonic behaviour in general. 

Such physical transonic accretion solutions can mathematically be reali zed as critical solutions on the phase portraits of the local ra- 
dial Mach number and the radial distance measured from the event horizon (Rav & Bhattachari eebool lA fshordi & Paczvnski 20031: 1511 
2003allbl:lRav & Bhattacharieel2005blla|Jchaudhurv et alhoOollRav & Bhattacharieel200dl2007al ; lBhattachariee & Ray|2007l ; lRav & Bhattacharied 
2007bl : lGoswami et alj| 2007; Bha ttachariee et alj|2009{j) 7 To maintain physical transonicty such critical points will perforce have to be sad- 



dle points, which will enable a solution to pass through themselves. In this connection, a "multi-critical" flow refers to the category of the 
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accretion configuration which can have more than one critical points accessible to the flow solution. For low angular momentum axisym- 
metric black hole accretion, it may so happen that the critical features are exhibite d more than once in the phase portrait of a stationary 



solution describing such flows, and accretion consequently becomes multi-critical (Liana & Thomson 1980; Abramowicz & Zurck 



B; irai et al]2004l : lAbraham et al.ll2 0P 
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In reality, such weakly rotating sub -Keplerian flows are in deed exhibited in various phy sical situations, such as detached binary systems 
fed b y accretion from OB stellar winds l llllarionov & Sunvaevl dl975h : lLiang & Nolanl (1984J)), semi-detached low-mass non-ma gnet i c bina - 
ries ( Bisikalo et al. 1 1998)), and super-massive black holes fed by accretion from slowly rotating central stellar clusters l llllarionov 1 1988); 
Holl ll999l) and references therein). Even for a st andard Keplerian accretion disc, turbulence may produce such low angular momentum flows 
(see, e.g., Igumenshchev & Abramowicz 1 1999b . and references therein). 



All of the aforementioned multi-critical flow dynamics are important in the astrophysical context. Such multi-critical behaviour allows 
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potentially provide an important and efficient mechanism for conversion of a significant amount of the gravitational energy into radiation by 
randomizing the directed infall motion of the accreting fluid. Shocks play an important role in governing the overall dynamical and radiative 
processes taking place in astrophysical fluid flows around black holes. 

Originally at a large distance, subsonic accretion encounters the outermost saddle type critical point and becomes supersonic. Subjected 
to the appropriate perturbative environment, such a supersonic flow encounters a shock and becomes subsonic again. The resulting flow 
has to pass through another saddle type critical point to meet the inner boundary condition as imposed by the event horizon. For accretion 
onto a black hole, the presence of at least two saddle type critical points is, therefore, a necessary (but not sufficient) condition for the 
shock formation. So multi-critical flow behaviour plays a crucial role in studying the physics of shock formation and related astrophysical 
phenomena. 

To understand the phase-space behaviour of low angular momentum shocked multi-transonic accretion, one usually constructs the 
corresponding autonomous dynamical systems analogue, and then identifies the saddle type critical points of the phase trajectories of the 
flow. Next, the global understanding of the flow topologies are performed which necessitates a complete numerical investigation of the 
non-linear stationary equations describing the velocity phase space behaviour of the flow. 

It is, however, still possible to semi-quantitatively realize the global behaviour of the transonic solution without taking resort to numerical 
techniques. Getting equipped with the mathematical formalisms of the general dynamical systems approach, it has recently been possible 
to conceive a clear ana l ytical conception of some of the global b ehaviours of the flow by analyzing the local features of the critical points 
1 Chaudhurv etap2006l : lMandal et al.ll2007l : lGoswami et alj2007h . 



It is important to note that along with understanding the critical point behaviour of the stationary accretion solution, it is also necessary 
to ensure that such stationary configuration are stable. This can be accomplished by studying the time evolution of a linear acoustic-like 
perturbation (around the stationary configuration) in the full time-dependent flow equations. Considering a hydrostatically balanced flow in 
vertical equilibrium, it has recently been observed that for accretion onto a non rotating black hole under the influence of various pseudo- 
Schwarzschild black hole potentials, t he characteristic feature s of the time development of the aforementioned perturbation ensures the 
stability of the stationary configuration dChaudhurv et al]|2006l) . The result obtained in this way was shown to be independent of the choice 
of the black hole potentials used to study accretion flow around a non rotating black hole. 

This work (Chaudhurv et al. I bood) . however, was performed for a particular type of flow geometry - hydrostatically balanced flow 



under the vertical equilibrium. Nevertheless, accretion processes onto astrophysical b lack holes are also studied for two othe r differen t flow 
geometries - flow with constant disc height, and flow under the conical eq uilibrium ( Liang & Thomson! 1980l : Abramowicz & Zureklll981 



Blaesll 19871 : iLu et al.lfl997l : Ichakrabarti & Dasll200ll : iGu & Foglizzol l2003l) (see section 2 for further detail about these two disc models). 
Those two flow geometries are relatively simpler to handle (in comparison to the flow configuration under the vertical equilibrium) without 
compromising the essential physics involved in the multi-transonic black hole accretion phenomena. In addition, these two flows are ap- 
propriate to study the low angular momentum inviscid flow configuration as well. It is thus instructive to investigate whether the stationary 
configurations remain stable (under perturbation) in these flow geometries as well. In other words, one needs to realize whether the stationary 
critical solutions are stable irrespective of the nature of the space time (choice of the black hole potential) as well as the flow geometry 
(structure of the accretion disc). 

This is precise objective of this work. The stationary and the time-dependent low angular momentum axisymmetric accretion around 
a Schwarzschild black hole, under the influence of a generalized pseudo-Newtonian black-hole potential in different flow geometries, have 
been analyzed. The stationary solutions have been considered to investigate their critical point behaviour, and to categorize systematically 
the nature of the critical points which appear in such flows. This is followed by a perturbative study of the full time-dependent flow, to follow 
the evolution of the perturbation and make predictions about the stability of the stationary configuration. Finally, observations have also been 
made about the nature of the acoustic metric embedded inside the flow. 
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2 THE EQUATIONS OF THE FLOW AND ITS FIXED POINTS 



When considering a rotating, axisymmetric, inviscid steady flow, the two most pertinent equations are the ones determining the drift in the 
radial direction (essentially Euler's equation), 



dv 1AP ,,. . A 2 „ 
v— + __ + 0'(r) - — = 
dr p dr r J 

and the equation of continuity, 



(1) 



(2) 



-(pvrH) = 0, 

in which, 4>(r) is the generalised pseudo-Newtonian potential driving the flow (with the prime denoting a spatial derivative), A is the conserved 
angular momentum of the flow, P is the pressure of the flowing gas and H = H(r) is the local thickness of the disc, respectively. The two 
foregoing equations give the steady continuum distribution of the velocity field, v(r), and the density field, p(r). But to close the two 
equations it will also be necessary to prescribe the functional dependences of both P and H on v and p, which, in the steady state regime, 
will imply an ultimate dependence on r. 

Following this requirement, the pressure, P, is first prescribed by an equation of state for the flow dChandrasekhan [l939). As a general 
polytropic it is given as P = Kp 1 , while for an isothermal flow the pressure is given by P = pnT/ (imn, in all of which, A" is a measure of 
the entropy in the flow, 7 is the polytropic exponent, k is Boltzmann's constant, T is the constant temperature, mu is the mass of a hydrogen 
atom and /i is the reduced mass, respectively. 

In fixing the function, H, one needs to look at the r elevant vertical geometry of the disc system. This can vary in many ways, with 



different degrees of complexity I Ch akrabarti & Das 2001). In the simplest case one could treat H to be just a constant, i.e. the disc is of 



uniform thickness. In the case of the conical flow dAbramowicz & Z urek 1981) one prescribes, H = Dr, where D is a simple constant 
of proportionality. While these two cases could be viewed as giving an explicit dependence of H on r, another well-invoked, but much 
more com plicated prescription in accretion literature is that of the disc with the condition of hydrostatic equilibrium imposed in the vertical 
direction 1 Matsumoto et alJl 1984| : Frank et alj 2002 ). In this particular instance, the function H in Eq. ® will be determined according to 
the way P has been prescribed 1 Frank et alj|2002 ). In all of these cases, however, it is a common practice to standardise transonicity in the 
flow by scaling its bulk velocity with the help of the local speed of sound, which is given as Cs = (dP/dp) 1/2 . With all these analytical 
requirements stipulated clearly, in what follows, the equilibrium properties of the flow will be studied for the three different kinds of vertical 
disc geometry mentioned above, under both polytropic and isothermal prescriptions for the equation of state. 



2.1 Polytropic flows 



With the polytropic relation specified for P, it is a straightforward exercise to set down in terms of the speed of sound, c s , the first integral of 
Eq. ([]} as, 



v 2 



r + 



2r 2 



£ 



(3) 



in which n = (7 — 1) 1 and the integration constant £ is the Bernoulli constant. Before moving on to find the first integral of Eq. {2]l 
it should be important to obtain the functional form of H. In simple cases of the vertical disc geometry, H usually becomes an explicitly 
defined function of the radial distance, H(r) — either constant disc height, or a conical profile with a linear dependence, H — Dr. Making 
a note of this fact, the first integral of Eq. lO could be obtained as 



4n 2 2 o-2 M 



(4) 



where M — (jK) n m 1 Abramowicz & Zurekll 1981 ) with m, an integration constant itself, being physically the matter flow rate. 

To obtain the critical points of the flow, it should be necessary first to differentiate both Eqs. (|3j and @, and then, on combining the two 
resulting expressions, to arrive at 



— (V ) = 

dr r 



\ 2 

— - r(p + c s 



H' 

l + r- 



(5) 



in which H' implies AH/Ar . The critical points of the flow will be given by the condition that the entire right hand side of Eq. ([5} will vanish 
along with the coefficient of d(u 2 )/dr. Explicitly written down, and following some rearrangement of terms, this will give the two critical 
point conditions as, 



r c <j>'(r c ) 



A 



1 + r c 



H'(r c 



H(r c ) 



(6) 



with the subscript c labelling critical point values. To fix the critical point coordinates, v c and r c , in terms of the system constants, one would 
have to make use of the conditions given by Eqs. ifSJ along with Eq. (O, to obtain 
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1 fj + i 

2 U-l 



r c 0'(r c ) 



1 + r c 



(7) 



from which it is easy to see that solutions of r c may be obtained in terms of A and £ only, i.e. r c = /i (A, £). Alternatively, r c could be fixed 
in terms of A and M . By making use of the critical point conditions in Eq. (4j one could write 



An 2 r 2 H 2 (r c ) 



r c <j>'(r c 



A 



1 + r c 



M 2 , 



(8) 



with the obvious implication being that the dependence of r c will be given as r c = /2(A, M). Comparing these two alternative means of 
fixing r c , the n ext logical step wo uld be to say that for the fixed points, and for the solutions passing through them, it should suffice to specify 
either £ or M Jchakrabartilll99(]|) . 

For the two relatively simple cases of disc geometry, i.e. constant H, and H — Dr, what Eq. (7J delivers are, 



1 fj+l 



2 \l - 1 
and 

' 7 + I 
4 U-l 



r c 0'(r c ) 



r c 0'(r c ) 



+ 4>(r c ) + 



+ (f>(r c ) + 



2ri 



2ri 



£. 



(9) 



(10) 



respectively. 

All the results obtained so far, can be compared with the case in which the vertical disc geometry is determined by the condition of 
hydrostatic equilibrium in the vertical direction. This requirement will deliver the functional form of H as, 

1/2 



H 



(11) 



which evidently shows that H is no more an explicit function of r. Now with the help of this form of H, and following the mathematical 
procedure outlined so far, the critical point coordinates could be fixed in terms of the system parameters, £ and M 2 , by the relations, 



27 
7-1 

and 

An 2 p 2 rl 



A 2 

r c <p'(r c ) - -j 



3-r, 



4>"{r c ) 



b'(r c ) P 2 



A 2 

r c <A'(r c ) - — 



3-r, 



4>"{r c ) 



4>>{r c ) 



2(n+l) 



(12) 



(13) 



respectively. A detailed presentation of all these re sults (pertaining to the sp ecific case of the disc balanced by hydrostatic equilibrium in the 
vertical direction) is to be found in a work done by Chaudhurv et al] 1 20061) . 



2.2 Isothermal flows 

For isothermal flows, the full mathematical treatment is actually much simpler. Here one has to go back to Eq. (TJ and use the linear 
dependence between P and p as the appropriate equation of state. On doing so, the first integral of Eq. (TJ is given as 

2 % 2 

— + c 2 lnp+^(r) + ^ =C (14) 

with C being a constant of integration. For flow solutions which specifically decay out to zero at very large distances, the constant C can be 
determined in terms of the "ambient conditions" as C = c 2 lnpoo. For the disc of constant thickness, or the disc with a conical flow, one is 
easily able to obtain an expression that is identical to Eq. (5j, from which the critical point conditions emerge exactly in the same form as 
what has been shown in Eq. (6j. However, major point of difference lies in the fact that in an isothermal system, the speed of sound, c s , is 
globally a constant, and so having arrive d at the critical point co nditions, it should be easy to see that r c and v c have already been fixed in 
terms of a global constant of the system 1 Chaudhurv et al.|[2006h . The speed of sound can further be written in terms of the temperature of 



the system as c s = ST 1//2 , where Q = (c/jimH) 1 ' 2 , and, therefore, it should be entirely possible to give a functional dependence for r c , as 
r c = /3(A, T). For the two cases of H being a constant, and H — Dr, once again the two respective relations for r c are 

A 2 

r c <t>'(r c ) J = c s ' ( 15 ) 

and 

A 2 

ro4>'(r c ) - — = 2c s 2 , (16) 
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with the subscript c labelling critical point values, as usual. 

To obtain similar results for the disc in vertical hydrostatic equilibrium, it is necessary to go back to Eq. dl It . with the restriction that 
7 = 1, and c s is a global constant. Then it becomes a simple exercise to fix the critical point coordinates, by deriving the expression, 



r c (t>'(r c ) 



A 2 



<t>"{r c ) 



<t>'(r.) 



(17) 



A detailed derivation of this particular result can once again be found in the work of Chaudhurv e t al. ( 2006). 



3 NATURE OF THE FIXED POINTS : A DYNAMICAL SYSTEMS STUDY 



The equations governing the flow in an accreting system are in general first-order non-linear differential equations. There is no standard 
prescription for a rigorous mathematical analysis of these equations. Therefore, for any understanding of the behaviour of the flow solutions, 
a numerical integration is in most cases the only recourse. On the other hand, an al ternative approach cou ld be made to this question, if 
the governing equations are set up to form a standa rd first-order dynamical system Jjordan & Smith! 19991) . This is a very usual practice 
in general fluid dynamical studies JBohr et al.ll 19930 . and short of carrying out any numerical integration, this approach allows for gaining 
physical insight into the behaviour of the flows to a surprising extent. As a first step towards this end, for the stationary polytropic flow, 
as given by Eq. Q, it sh ould be necessary to parametrise this equation and set up a coupled autonomous first-order dynamical system 



. given by bq. it sh 
, (Jordan & Smith! 19990 



2v 2 



A 2 ±l 2 / H' 
— - r4> + Cs I 1 + r— 



dr 

d7 



= r I v 



(v 2 - Cs) 



(18) 



in which r is a n arbitrary mathematical parameter. With respect to accretion studies in particular, this kind of parametrisation has been 
reported before feav & Bhattachariedl2002l ; kfshordi & PaczyfiskilliooH ; Ichaudhurv et alJbood : iMandal et ai]l2007l ; iGoswami et~ail2007 ; 
Bhattachariee et al.ll2009blv 



The critical points have been fixed in terms of the flow constants. About these fixed point values, upon using a perturbation prescription 



of the kind v 



■ + Sv 2 



c sc + Sc s and r = r c + Sr, one could derive a set of two autonomous first-order linear differential equations 



in the Sr — Sv 2 plane, with 5c 2 having to be first expressed in terms of Sr and Sv 2 , with the help of Eq. ([4} — the continuity equation — as 



5c 2 , J^ 2 



l + r c 



H'jvc) 
H(r c ) 



(19) 



The resulting coupled set of linear equations in Sr and Sv' 2 will follow simply as 

H'(r c )' 
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H'(r c 



H(r c ) 



Sr . 



(20) 



in which a prime implies a derivative with respect to r. Trying solutions of the kind Sv 2 
deliver the eigenvalues fl — growth rates of Sv 2 and Sr — as 

2 



exp(Slr) and Sr ~ exp(ilr) in Eqs. ( 120b . will 
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H(r. 
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+ <t>'(r c ) + r c <t>"(r c ) + (7 - 1) 1 1 + r 



H'(r c ) 
' H(r c ) 
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For the specific cases of the simple vertical geometries, i.e. H is a constant, and H = Dr, the foregoing expression reduces to 

v2 „2 



n =(7-1) ct, 

and 



(7 + 1) r c ct 



^4r + 4>'(r c ) + r c 0"(r c ) + (7 - 1) ^ 
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(23) 
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Figure 1. Different regions in parameter space of E and A. corresponding to number and nature of critical points for polytropic flows. The dotted lines are for 
vertical equilibrium geometry, the solid lines are for conical geometry, and the small dashed lines are for constant-height disc. For notations see text. 



respectively. 

A similar treatment might also be extended to the case of the disc in vertical hydrostatic equilibrium. Following a perturbative treatment 
about the fixed point coordinates, the eigenvalues, f2, could be obtained in this instance as 
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(24) 



where Af^r) = r 3 (f>' (r), and with 



A = r, 



4>'{r c ) 



B = 1 + r, 



<l>"'(r c ) _ 4>"(r c ) 
<t>"{r c ) C <t>'{r c ) 



The detailed calculations to arrive at Eq. d24b could once again be accessed in the work of lChaudhurv et al. I (hoodi . 

For isothermal flows, similar expressions for the related eigenvalues may likewise be derived, given a particular form of the function, 
H(r). The algebra in this case is much simpler and it is an easy exercise to assure oneself that for isothermal flows one just needs to set 7 = 1 
in Eqs. J21b . J22I J23b and d24b . to arrive at a corresponding relation for fl 2 . However, it should be incorrect to assume that in this kind of 
study, one could always treat isothermal flows simply as a special physical case of general polytopic flows. For polytropic flows, the position 
of the fixed points, under a given form of 4>(r), will be determined by Eq. or by Eq. (8), in the case of simple vertical disc geometries, 
and by Eq. l !12b or by Eq. J13b , in the case of the disc balanced by hydrostatic equilibrium in the vertical direction. On the other hand, for 
isothermal flows the fixed points are simply to be determined from the critical point conditions themselves, since c s is globally a constant in 
this case 1 Chaudhurv et al The resulting difference is by no means trivial. 

Once the position of a critical point, r c , has been ascertained, it is then a straightforward task to find the nature of that critical point by 
using r c in either Eq. ( 121b or Eq. d24b . depending on the disc geometry. Since it has been discussed in Section[2]that r c is a function of A and 
T for isothermal flows, and a function of A and £ (or M) for polytropic flows, it effectively implies that Q 2 can, in principle, be rendered as 
a function of the flow parameters for either kind of flow. A generic conclusion that can be drawn about the critical points from the form of 
Q 2 in Eqs. ( 121b and ( 124b . is that the only admissible critical points will be saddle points and centre-type points. For a saddle point, Q 2 > 0, 
while for a centre-type point, Q 2 < 0. 



4 NUMERICAL RESULTS 

It is quite evident from previous discussions that the location of critical points and their nature (saddle or centre) are easily determined from 
the roots of the equations (Eq 191 1111511 7b and the sign of Q 2 . These, in turn, effectively extract the qualitative features of the phase portrait 
without having an explicit plot of it. 

In this work, with numerical evidence, these tools are exploited to compare the features of the accretion flow in different disc geometries 
under polytropic as well as isothermal equations of state. Though this fo rmulation of the procedure works for any choice of a pseudo- 
Schwarchild potential, the Paczynski-Witta potential 1 Paczvhski & J.fl98(il) . 



<f>pw(r) = - 2(r 1 _ ^ , (25) 

has been selected for the presentation of numerical computations within pseudo-Newtonian framework. 

First of all, Fig. Q] shows the variation of the number of critical points for polytropic flows with the variation of £ and A for the three 
disc geometries. For example, the region bounded by Ah, Bh and Dh depicts the parameter values corresponding to three critical points in 
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Figure 2. Different regions in parameter space of Tin and A corresponding to the number and nature of critical points for polytropic flows. Tio means 
temperature of the flow in 10 10 °K. Linestyles are similar to the previous figure. For other notations see text. 
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Figure 3. Variation of Q 2 with A at constant E for polytropic flow. 



the constant energy flow solutions, for the constant-height disc geometry, within which for the subregion bounded by Ah, Bh, and Ch the 
value of M at the innermost critical point (Mi n henceforth) is higher than that at the outermost critical point {Mouti- The relative values of 
M are just the reverse for the other subregion bounded by Ah, Ch and Dh- On the curve AhCh at all points A4 acquires the same value 
at both the critical points; for all other regions shown in the graph there will exist only single critical points. 

For other disc models, similar features are depicted by the symbols with suffixes C (Conical) and V (Vertical equilibrium). Hence for 
all three geometries, there are certain wedge shaped regions which correspond to three critical points and outside these regions, the parameter 
values within the parameter space shown in the figure, generate single critical points only. Hence for low values of £ with there will be an 
interval of values of A for which three critical points exist (opening up the the possibility of multi-transonic accretion within the subregion 
where A4 in > A4 out under certain other condition (see later)). In Fig.[2]the same thing is shown with the relevant parameters T and A. Here 
also the appearance of the wedge shaped region conforms to the possibility of multi-transonic accretion for certain interval of values of A. 
But here in the c 2 — Q 2 T — > limit, all three equations merge together and as a result of it, in the figure the lower points B and C for 
all models, merge together. Instead of M, here C is the quantity which makes a difference between subregions ABC and ACD in different 
models. 

Both in polytropic and isothermal cases, the interval of A permitting more than one critical points progressively shrinks for higher values 
of the other parameter (i.e. £ and T, respectively) and there are certain values of £ or T, dependent on the disc geometry, over which there 
will be only a single critical point, whatever be the values of A. The wedge shaped regions also widely vary from one disc geometry to 
another. Actually the regions of parametric values allowing multicritical points for various disk geometries may be approximately mapped 
from one another with some suitable scaling of 7 (( Chakrabarti & Dasl2001 )); but if one may map such a parametric region for the polytropic 
flow into another correspondoing parametric region with a choice of 7 (for the second case) nearly equal to unity, it should not create the 
impression that the first flow in the polytropic condition may at least approximately mimic the second flow in the isothermal condition, 
because the physical conditions and the first integral of motions under the two conditions (polytropic and isothermal) are distinctly different. 

The nature of the critical points will be apparent from the Fig.[3]and Fig. [4] In these plots a positive value of Q, 2 indicates a saddle point 
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Figure 4. Variation of Q with A at constant T for isothermal flow. 



while a negative value does the same for a centre. A common feature of both the separate plots is that initially there is region of a single saddle 
point, followed by the birth of a centre-type point and another saddle point (saddle-centre type bifurcation); then at another higher value of 
A, the centre-type point coalesces with the other saddle point (the outer one) and both of them annihilate each other (another saddle-centre 
type bifurcation, but this time in the opposite direction) so that the remaining saddle point (outer) "survives" above the critical value of A. 
Actually in the underlying phase plots (not shown here), among the three critical points, there will always be a pair of centre-saddle for which 
the separatrices of the saddle form a homoclinic connection around the centre, as long as the other constant of flow {M for polytropic case 
and C for isothermal case) discriminates between the two saddle points. The remaining saddle point will have the separatrices connecting the 
event horizon with the infinity. The separatrices of the saddle point with higher value of M in the polytropic case (with lower value of C in 
isothermal case) form the homoclinic connection. Physically this is what it has to be, because it is the third critic al point (saddle) wh i ch wi ll 



allow the transonic flow solution from infinity to the event horizon, and here, following the line of argument in Rav & Bhattachariee (2002), 
it can be stated that the stationary flow solution has to settle on the separatices of a saddle point only. That is why among the critical points, 
only the saddle points are termed as sonic points, and not the centre-type points, although for both the types of critical points the flow speed 
is equal to the sound speed, i.e. Mach number becomes unity. 

So if the value of A is increased from a sufficiently low value, first there will a span of single saddle point (solid line in the figure) 
then the bifurcation occurs resulting in a pair of centre-type points (with intermediate value r, denoted by r m id) and a new saddle point 
(location Tin) so that the older saddle point becomes the outermost critical point (location r out ). Before A attains a certain value, the relation 
Min > Mout (for isothermal case d„ < C ou t) is maintained and for this interval of A, the Q 2 functions for the three critical points are 
plotted with dotted lines. Here the separatrices of the newly formed saddle point encompasses the centre-type point, forming a homoclinic 
connection; after exceeding this value of A, what happens is that Min > Mout, and the homoclinic connection for the inner saddle opens 
up, giving rise to the same sort of homoclinic connection from the outer saddle. Within this phase, the function of Q, 2 is depicted with small 
dashed curves. At the end of this domain of A, the centre-type and the outer saddle point coalesce and "destroy" each other, after which, only 
the inner saddle point exists with increasing value of A (solid curves). 

The scenario in the isothermal case is almost similar. There is an interesting situation along the boundary between the two subregions 
within the wedge shaped region corresponding to multicritical points. For these parameter values Min = Mout, and thus two saddle points 
will be linked by a heteroclinic connection of their separatrices; for the A just below this value there is some homoclinic connection for the 
inner saddle point, and just above this value the homoclinic connection corresponds to the outer critical points. From the point of view of a 
general dynamical system along this boundary, a type of bifurcation occurs that may be termed as a heteroclinic one. 

The boundaries of the region in the parameter space permitting multicritical points, as it has been discussed already, are associated with 
saddle-centre bifurcation or merging of a pair of roots of the equations (Eg 1911 ll I15ll7t needed for determining the critical points under 
various conditions discussed above. Now all these equations are polynomial equations. The discriminant of a general polynomial, 

P n (x) = a n x" + a„_ii™ _1 H + aix + a , (26) 

can be expressed as in terms of its roots, Xi's, as 



Xi-Xj) 2 . (27) 



The discriminant may be expressed as the determinant of a matrix called Sylvester matrix (http://mathworld.wolfram.com/PolynomialDiscriminant.html 
and other references therein), 
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(28) 



up to a factor. Hence it is expected that det (S) will be zero on the abovementioned boundaries and actually it is so. Here the contour plot of 
det (S) — value for the polytropic flow in constant height geometry (i.e. for the polynomial in r c in Eq[9j in S-X space is shown in Fig|5] 
The curve exactly conforms with the corresponding boundary curve (in small dashed style) in Fig[TJ drawn on the basis of the numerical 
method of detailed root finding. So this procedure may be thought of as a much easier alternative to find the multicritical parametric values. 

Considering now the question of actual accretion process, from Eq|5]it is evident that for conical and constant height model the critical 
points are actual sonic points, though these two are not exactly same for the other model.As for all parameter values there will always be one 
critical (sonic) point with separatrices spanning from infinity to the event horizon, there will always exist a transonic solution. But the features 
of the critical points in all the three models opens the possibility of multitransonic accretion too within the region with multicritical points 
under certain conditions. First, there should be some transition from the supersonic branch of the non-homoclinic separatrices (trajectory) 
of the critical point (sonic point) with global transonic solution to the subsonic branch of the homoclinic solution of the other sonic point at 
a particular value of r, Secondly there should be some quantity which will discriminate between the trajectories of these two sonic points. 
For such a transition one such possibility is the occurrence of a shock, which may take place obeying the condition Mm > Mout for the 
polytropic case (and Cm < C ou t for the isothermal condition). After making the transition, the flow will pass through the corresponding inner 
sonic point. The determination of the exact location and the exact parameter values for such a transition really taking place is out of the scope 
of the present paper, but here it is touched upon to indicated that this full analysis keeps alive the possibility for such types of multitransonic 
accretion, irrespective of the disc geometries and equations of state. 



5 TIME-DEPENDENT STABILITY ANALYSIS OF STATIONARY SOLUTIONS 

The time-dependent generalisation of the continuity condition for an axisymmetric pseudo-Schwarzschild disc is given as 

^ + i^(E«r) = 0, (29) 
at r or 

in which the surface density of the disc, S, is to be expressed as £ = pH dFrank et al.l l2002). From the foregoing expression, one can, 
therefore, obtain 



dp Id, ... 



(30) 



Defining a new variable / = pvrH, it is quite obvious from the form of Eq. l |29t that the stationary value of / will be a con stant, ,fo, which 
can be closely identified with the matter flux rate. This follows a similar approach to spherically symmetric flows made bv lPetterson et al 
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1 1980h and lTheuns & David 1 199 lh . For an axisymmetric di sc, which ha s no dependence on any angle variable, this approach has also been 
adopted for a flow driven simply by a Newtonian potential (Ray 2003c). The present treatment, of course, is of a more general nature, the 
disc flow being driven by a general pseudo-Newtonian potential, <f>(r) for all disk geometries. In this system, a perturbation prescription of 
the form v(r, t) = Vo(r) + v'(r, i) and p(r, t) = po(r) + p'(r, t), will give, on linearising in the primed quantities, 



fo Po v 



(31) 



which is a relation that connects all the three fluctuating quantities, v' , p and /', with one another (here the subscript denotes stationary 
background values in all the cases). Going back to Eq. J30t . it becomes possible to connect p' exclusively to /', through the relation 



dp' vqPq 
dt + f 



df 
i:)r 



0. 



(32) 



The cas e of the disc being balan ced vertically under hydrostatic equilibrium is very different in mathematical terms, and has been taken 
up earlier by Chaudhu rv et al. N2006h . Nevertheless it is worth going back to it for making some interesting comparisons. By using Eq. i ll It 
and the polytropic relation P = Kp 1 ', Eq. d29t can be rendered as 



d 
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from which, under a new definition, / = p ( - 7+1 ' //2 «r 3 ^ 2 /\ 
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From Eq. d33t , it is also very easy to set down the density fluctuations, p' , in terms of /', as 

3 2 VoPO 



dt 



fo 



dr ' ' 



(33) 



(34) 



(35) 



with f3 2 = 2(7 + 1)~ , as before. This result may be compared with Eq. d32b and the difference noted. If, however, one were to study 
an isothermal flow balanced by hydrostatic equilibrium in the vertical direction, then Eq. dl It would have to be constrained by 7 = 1 and 
c s being constant. Under these conditions, the expression for density fluctuations in the flow will be identical to Eq. d32b . rather than be 
described by Eq. d35t . 

The equation for density fluctuations may show variations in terms of a constant scaling factor under different vertical disc height 
geometries, but no such thing happens for the velocity fluctuations. Combining either Eqs. ( 131b and d32t . or combining Eqs. J34t and Q51 > 
will render the velocity fluctuations as 



dv' 
~dt 



df , df 

fo \ dt or 



which, upon a further partial differentiation with respect to time, will give 
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ffl] 



The time-dependent equation for the radial drift is given as 

A 2 

from which the linearised fluctuating part could be extracted as 



dv ^ dv ^ 1 dP _|_,//% 
dt dr p dr 
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(36) 



(37) 



(38) 



(39) 



with Cso being the speed of sound in the steady state. Differentiating Eq. l |39t partially with respect to t, and making use of either Eq. l |32t or 
Eq. l |35t , along with Eqs. d36l l and l |37| l, to substitute for all the first and second-order derivatives of v' and p', will deliver the result 
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(40) 



in which either a = 1 or a — f3 2 , depending on the choice of a particular disc geometry and the equation state applied. For isothermal flows, 
(7 = 1, for whatever disc geometry one considers — H is constant or H — Dr or H is as it is described by Eq. QjJ , The same value of a is 
also obtained for polytropic flows in the first two simple cases of the height function, H. The common feature running through all these cases 
is that H in Eq. d30t does not have any dependence on time. It is only when the flow is polytropic and the disc height geometry is expressed 
by Eq. dl It . will H have a time-dependence, whose ultimate consequence will be that a = /3 2 in Eq. d40t . 
All the terms in Eq. d40t can be expediently rendered into a compact formulation that looks like 
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(41) 



in which the Greek indices are made to run from to 1, with the identification that stands for t, and 1 stands for r. An inspection of the 
terms in the left hand side of Eq. d4QI > will then allow for constructing the symmetric matrix 



f" 



Vo 



fo \v Vq- G(? sQ 



Now the d' Alembertian for a scalar in curved space is given in terms of the metric g M „ by l Visserl 19981) 

with g M " being the inverse of the matrix implied by gp, u . Using the equivalence that f M " = ^/— gg M 
immediately possible to set down an effective metric for the propagation of an acoustic disturbance as 



(42) 



(43) 



Soff 



Vq 



Vo 



and therefore g = det (P"), it is 



(44) 



which can be shown to be ent irely identica l to the metric of a wave equation for a scalar field in curved space-time, obtained through a 
somewhat different approach l Vissei 1998h . The inverse effective metric, gj^, can be easily derived by inversio n of the matr ix given in 
Eq. l !44t . and this will give v 2 = <7Cg as the horizon condition of an acoustic black hole for inflow solutions ( IVisserlll998r) . From the 
perspective of the propagation of acoustic waves (carrying information in any fluid system) in the accretion disc, what can be concluded from 
these arguments is that the disc geometry and the equation of state act together in determining the speed of information propagation. In the 
case of the disc being supported by hydrostatic equilibrium in the vertical direction, if the flow is polytropic, then the speed of information 
propagation will be less than the speed of sound by a factor, /3. So transonicity will not take place exactly when the bulk flow speed, v, 
becomes equal to the speed of sound, c s . In all the other cases, transonicity will, however, take place when v = c s . 
Finally, a little readjustment of terms in Eq. d40t will give an equation for the perturbation as 
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(45) 



whose detailed solution has been given in earlier works on inviscid axisymmetric flows dRavH2003c ; Chaudhurv et allipod) 



6 CONCLUDING REMARKS 

The primary motivation behind this work, as has already been mentioned, is to study the transonic properties and the stability issues of 
low angular momentum axisymmetric black hole accretion from a dynamical systems point of view. Flows around a non rotating black 
hole in three possible geometries, namely, the constant height, the conical equilibrium and the hydrostatically balanced vertical equilibrium 
respectively, have been studied under the influence of a generalized pseudo-Schwarzschild black hole potential. Both polytropic as well as the 
isothermal equation of state have been used to describe the flow. The analytical formalism provided here ensures that even without performing 
explicit numerical integration along the streamlines, a reasonably clear picture of the behaviour of the phase trajectories of the flow can be 
well apprehended. This work, thus, can be considered as a crucial alternative approach to understand the behaviour of non-exactly solvable 
coupled differential equations describing autonomous dynamical systems, as they are to usually encountered while studying the accretion 
processes around compact astrophysical objects. 

The stationary critical solutions, as well as the stability criteria of such stationary configurations have been analyzed. For stationary 
structure, one finds that the general profile of the parameter space marking the multi-critical flow solution remains unaltered for all three flow 
geometries, as well as for every black hole potential available in the literature. The corresponding numerical domain of [£, A, 7], for which 
the multi-criticalit y has been observed , is, however, different for different disc geometries, as well as for different potentials for a particular 



equation of state. Chakra barti & Dasl a200lh proposed that the three different disc models, namely constant height flow, flow in conical 
equilibrium and in vertical equilibrium should be identical, provided one can change the polytropic constant in certain ways. In this way 
they suggested that a relativistic flow of constant height might have same properties as an isothermal flow in two other geometries. This is, 
however, not true. One is able to argue that it is not possible to allow accretion under various flow geometries identically by merely changing 
the polytropic index of the flow, since the shift from one flow geometry to another is a non-continuous transformation in a multi-dimensional 
parameter space. Hence, it is certainly not a mathematically consistent approach to establish an equivalence between the polytropic flow in 
one type of disc geometry to the isothermal flow in the other by making a 7 transformation, since the corresponding first integrals of motion 
for the flow, derived on using the two different equations of state (polytropic/adiabatic and isothermal) imply distinctively different physical 
properties. 

The Mach number relation at the critical point differentiate the flow in vertical equilibrium (M c < 1 for polytropic flow) from the two 
other kinds of flow (constant height flow and the flow under conical equilibrium), for which the critical point and the sonic point are identical 
(i.e., M c = 1), upon using both the equations of state. The result M c < 1 is a consequence of the effective speed of the propagation of the 
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linear acoustic perturbation (sound speed) due to the presence of the polytropic sound speed c s in the expression of the disc height (which, 
in other words, is a direct repercussion of the vertical equilibrium assumption itself). 

The stability properties of the aforesaid stationary configuration have been realized by perturbing (about the stationary configuration) the 
full time-dependent flow equations in various disc geometries, under various black hole potentials, and then by observing the time evolution 
of such perturbations. The general form of the wave equations corresponding to the dynamics of such perturbations, as well as the related 
acoustic metric, are identical for any black hole potential used, but is very different for different flow geometries. However, such perturbations 
do not diverge in any physical sense for any kind of flow geometries. This ensures that the stability of the stationary configuration, at least 
for accretion around non-rotating black holes under various pseudo-Schwarzschild potentials, is a generic feature independent of the nature 
of the space time (the explicit form of the black hole potential) as well as the geometric configuration of the flow (disc structure). 

The exact forms of the corresponding acoustic metric for various disc geometries have also been derived in this work. This will help 
in studying the low angular momentum pseudo-Schwarzschild axisymmetric black hole accretion as a natural example of analogue gravity 
phenomena. It has recently been shown that the a nalogue surface gravity can be computed for multi-critical accretion onto astrophysical 
black holes i Abraham et al. 20061 : Das et al. 2007 ). The present work will allow to study such phenomena for various different geometric 
configuration of the flow. 
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